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The recently announced Higgs discovery marks the dawn of the direct probing of the electroweak 
symmetry breaking sector. Sorting out the dynamics responsible for electroweak symmetry breaking 
now requires probing the Higgs interactions and searching for additional states connected to this 
sector. In this work we analyze the constraints on Higgs couplings to the standard model gauge 
bosons using the available data from Tevatron and LHC. We work in a model-independent framework 
expressing the departure of the Higgs couplings to gauge bosons by dimension-six operators. This 
allows for independent modifications of its couplings to gluons, photons and weak gauge bosons 
while still preserving the Standard Model (SM) gauge invariance. Our results indicate that best 
overall agreement with data is obtained if the cross section of Higgs production via gluon fusion is 
suppressed with respect to its SM value and the Higgs branching ratio into two photons is enhanced, 
while keeping the production and decays associated to couplings to weak gauge bosons close to their 
SM prediction. 



I. INTRODUCTION 



The electroweak symmetry breaking (EWSB) mechanism has been elusive for many decades. However the recently 
announced discovery of a 125 GeV Higgs boson at the CERN Large Hadron Collider (LHC) 0, Hi opens a new 
era in particle physics. The pressing questions now are related to the properties of this new observed state, like its 
spin and couplings, in order to extend our knowledge of the EWSB sector. In this work we employ the data used for 
the Higgs discovery to constrain its couplings to gauge bosons. 

Presently there are many possible EWSB scenarios ranging from the Higgs being elementary and weakly interact- 
ing , as in the Standard Model, to it being composite and related to a new strongly interacting sector [HI . In 
this last case the precision electroweak measurements and flavor changing neutral currents lead to strong constraints. 
However, recent theoretical advances have made possible the construction of models in agreement with the experi- 
mental bounds [T2j . The distinction between the different scenarios can be carried out by looking for further new 
states associated with the EWSB mechanism and/or by careful studies of the Higgs boson couplings. 

In this work we assume that the observed Higgs boson is part of a SU (2) l doublet and that possible additional states 
are heavy enough not to play a direct role in the low energy phenomenology. This assumption is realized in models 
where the Higgs boson is a pseudo-Goldstone boson of a larger broken global symmetry |13l - [l9f . Under this assumption 
we consider the most general dimension-six effective Lagrangian invariant under linear SU{i)c ® SU{2)l (8i U{1)y 
transformations to describe the interactions of the Higgs boson with the electroweak gauge bosons, as well as with 
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the gluons [20, l2l| . For the sake of simplicity we assume that the Higgs has the same interaction with fermions as in 
the SM, nevertheless this hypothesis still has to be tested further^. This scenario can be falsified by the discovery of 
new states or by the non-observation of its predictions to the triple electroweak-gauge-boson vertices. 

The effective operators describing the Higgs anomalous interactions modify both the Higgs production mechanisms 
and its decay patterns, therefore we combine several channels to unravel the contribution of the different operators. In 
our analyses we use the most recent data from the Tevatron [HI and LHC at 7 TeV [Ullll and at 8 TeV [IHlli^. 
Anomalous interactions also enhance the Higgs decay into as well as its production in association with a photon. 
Nevertheless, the available statistics is not enough to make these channels visible. 

This article is organized as follows. In Section II we introduce the dimension-six effective operators and the different 
scenarios studied in this work. Details of our analyses are presented in Section HI and Section IV contains their results. 
Finally we discuss the main conclusions in Section V. 



II. HIGGS ANOMALOUS INTERACTIONS 



In this work we assume that even if there is new physics associated with the electroweak symmetry breaking sector, 
the Higgs boson observed at LHC is still part of a SU{2)l doublet, the SM gauge invariance holds and no additional 
light states, relevant to the Higgs observables, are present in the spectrum. Under these assumptions the new effects 
can be parametrized in a model independent way by extending the SM with the addition of higher dimension operators 
that are invariant under linear SU{3)c <8> SU{2)l U{1)y transformations. 

In this framework the first corrections to the Higgs couplings to gauge bosons are expressed as dimension-six 
operators that can be written as 

= E ' (1) 

n 

where the operators 0„ involve vector-boson and/or Higgs-boson fields with couplings /„ and where A is a charac- 
teristic scale. Requiring the operators 0„ to be P and C even, there are only seven dimension-six operators that 
modify the Higgs-boson couplings to electroweak vector bosons and one to gluons [1^ [2l| : 

Ogg = G^l.G'^^'' , Oww = $^W'p.VF'"'$ , Obb = $+5^.5"-$ , 

Obw = <^>^B^,W^"'^ , Ow = (i^^$)tW^^'^(i^,$) , Ob - (I?^$)t-B^'^(i?,$) , (2) 

= p^$)^$t$pM$) , 0$^2 = (^t*) 9^ ($t$) , 

where $ stands for the Higgs doublet, is the covariant derivative, 13^^^ = i{g' /2)B^^y and W^j,y = i(g/2)cr°W°j,, 
with B^i,, W^^, and G"^, being respectively the [/(l)^, SU{2)l and SU{3)c field strength tensors. We denote the 
SU{2)l {U{1)y) gauge coupling as g (g') and the Pauli matrices as a"^. 

The effective operators in Eq. ^ give rise to anomalous Hgg, Hjj, HZj, HZZ, and HW'^W~ couplings, which 
in the unitary gauge are given by 

= 9Haa HGl^G'^'^'' + gn^, HA^^A'^'^ + g^^\^ A^,Z>^d"H + g^^\^ HA^^Z^^'' 

+ 9^HZZ Z^^a^d^B + .gL'L ^^p.^''" + 9h\z HZ^^Z^^ (3) 
+ .9ww {W+'^W- ^d^B + h.c.) + g%^ BW+W' + g%^ BW+W- , 

fl 2l fl 2 3) ('12 3) 

where V^^ = - d^Vfj, with V = A, Z and W. The effective couplings gngg, gn^yy, gnzj' 9hww ^^'^ 9hzz are 



^ The preliminary CMS results indicate that the SM values for the Higgs couplings to fermions are within the 90-95% CL allowed 
region. 
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related to the coefficients of the operators appearing in ([Ij through, 

foGV 



9Hgg 



yHZi 



(2) 



9hzz 



(2) 

9hzz 



(3) 

9hzz 



9hww 



(2) 

9hww 



(3) 

9hww 



as fgV 

A2 ~ Stt A2 ' 
/ gMw\ s^ifBB + fww ~ Ibw) 



V A2 y 2 

gMw\ s{fw - fs) 



A2 J 2c 

gMw \ s[2s^fBB - 2c^fww + (c' - s^)fBw] 
A2 J 2c 

gMw\ c^fw + s^fs 



A2 ; 2c2 

gMw\ s'^fsB + c^fww + c^s^Jbw 



A2 ; 2c2 
gMwv'^\ Us - h.2 



(4) 



A2 ; 4c2 

9^Jw\ fw 
A2 J 2 ' 

9]\Iw\ , 

I fww , 

A2 



where s and c stand for the sine and cosine of the weak mixing angle respectively. We notice that we have rescaled the 
coefficient Jgg of the gluon-gluon operator in terms of a coupling fg also including a loop suppression factor. In this 
way an anomalous gluon-gluon coupling fg ^ 0(1 — 10) gives a contribution comparable to the SM top loop. For the 
operators involving electroweak gauge bosons we have kept the normalization commonly used in the pre-LHC studies, 
for example, in Refs. [28l - [3^ . The couplings g'^^zz ^'^'i 9'hww include the effects arising from the contribution of the 
operators 0$,! and C'$,2 to the renormalization of the weak boson masses and the Higgs field wave function. 

For the sake of concreteness in this work we focus our attention on modifications of the Higgs couplings to gauge 
bosons associated with the five operators Oggi Obb, Oww, Ob, and Ow- The operator Obw contributes at tree 
level to the W^-B mixing and is therefore very strongly constrained by the electroweak precision data [28l l29l l33l |34{ . 
Similarly 0$.i contributes to the Z mass but not to the W mass and it is severely constrained by the p parameter. 
Moreover the operators O^^ ^'^d ©$,2 lead to a multiplicative contribution to the SM Higgs couplings to ZZ and 
WW . Thus in the present analysis we do not consider effects associated with Obw, 0$,! and 0^^2 as their coefficients 
are already very constrained or their possible effect on the measured Higgs observables is degenerated with that of 
the five operators considered. Their impact on the Higgs phenomenology can be seen in Refs. [sl. |35| - |40| . 

Notice also that one expects the contribution of new physics to the five operators considered to take place at loop 
level [4l| . Therefore, we expect that the largest effect of these effective interactions should appear in the couplings of 
the Higgs to photon-photon and gluon-gluon since these couplings take place through loop effects in the SM. 

One important property of the operators Ob and Ow is that they also modify the triple gauge-boson couplings 
^W^W~ and ZW~^W~. Consequently they can be directly probed in additional channels not directly involving the 
Higgs boson [3ll l43j . The triple gauge-boson effective interaction can be rewritten in the standard parametrization 
of the C and P even interactions 143 : 



-wwv 



-i9wwvi^9Y 



w+v,w- 



KV 



w+w-v" 



-w+ w- 



(5) 



where gww-y = e and gwwz = e/{sc). In general these vertices involve six dimensionless couplings gY , kv, and Ay 
(y = 7 or Z). Notwithstanding the electromagnetic gauge invariance requires that gj = 1, while the remaining five 
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free couplings are related to the dimension-six operators that we are considering: 





J 1 


1 = 


I ml 

IW , 




— — hi-y ' 


1 = 


2 A2 l^'^ ^ 


Akz 




1 = 


1 ™| / 2 . 




= Xz 




. 



(6) 



In summary, in the theoretical framework that we are using the observables depend upon 5 parameters, /g, fs, fw, 
Jbb and fww- In what follows for the sake of simplicity we focus on two different scenarios: 

• Scenario I: we impose that fw = fs and fsB = fww- This scenario has three free parameters {fwi fww 
and fg) and it exhibits a constraint between the three couplings of the Higgs to electroweak vector bosons. This 
scenario predicts the existence of anomalous triple electroweak gauge-boson interactions. 

• Scenario II: we set fw — fs — and fww = fsB- This scenario has two free parameters {fg and fsB — fww) 
and it can be considered the low-energy limit of an extension of the SM that contains an extra heavy scalar 
multiplet; for details see Ref. (s^). Moreover, this scenario cannot be constrained by data on triple gauge-boson 
couplings. 

The above relations ([6]) allow us to constrain the couplings fs and fw using the available experimental bounds 
on the effective couplings Agf , Akj and Akz [11]. Nevertheless these experimental bounds are usually obtained 
assuming only one anomalous operator different from the SM value at a time, an assumption which is not consistent 
with our scenario I. For this reason strictly speaking one cannot apply the exclusion limits in Ref. (45j to this scenario. 
Nevertheless if we assume no strong cancellations between the contributions of the different triple gauge-effective 
operators we can estimate the size of the exclusion limits on fw and fs- Using the 95% CL regions from Ref. [45| on 
Agf , Akj or Akz we obtain that in scenario I the 95% CL regions on fw = fs are [—13, 7] TeV"^, [—18, 9] TeV"^, 
and [—85, 20] TeV"^ respectively. Notice also that, LHC already with present runs has potential to constraint the 
triple-gauge boson vertices [i^ and the collaborations are starting to look for deviations [12,111]. However, at present, 
their individual limits have not reached the level of the LEP bounds yet. 



III. ANALYSES FRAMEWORK 



In order to obtain the present constraints on the Higgs anomalous interactions we perform a chi-square test using 
the available data on the signal strength (/x) from Tevatron, LHC at 7 TeV and LHC at 8 TeV. We assume that the 
correlations between the different channels are negligible except for the theoretical uncertainties which are treated 
with the pull method [U 113 in order to account for their correlations. 

Schematically we can write 



X = mm > ^ h > — — 7) 

J ■> pull ' 



where j stands for channels presented in Tables U and HIl We denote the theoretically expected signal as /ij , the 
observed best fit values as fi^^^ and errors as f^'"- As we can see from these tables the errors are not symmetric, 
showing a deviation from a Gaussian behavior as expected from the still low statistics. In our calculations we make 
the errors in each channel symmetric by taking 



Concerning the theoretical uncertainties, the largest are associated with the gluon fusion subprocess and to account 
for these errors we introduce two pull factors, one for the Tevatron (^y) and one for the LHC at 7 TeV and LHC 
at 8 TeV (^l). They modify the corresponding predictions as shown in Eqs. ([T2|) and ((T3|). We consider that the 
errors associated with the pulls are ar = 0.4 and (7^ = 0.15. As statistics build up it will be necessary to introduce 
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Channel 


exp 


Comment 






CDF & D0 [22] 


pp — >■ bb 


2.01^7 


CDF & D0 [22] 


pp — >■ 77 


3.6t^;° 


CDF & D0 [22] 


pp TT 




ATLAS [23] 


pp — >■ bb 


U.0_2 


ATLAS ^ 


pp zz* e+re+r 




ATLAS [23] 


pp^ww* ^i+vrv 




ATLAS [23] 


pp — >■ 77 


r, r, + 0.8 

^•^-0.8 


ATLAS [27] 


pp — >■ TT 


0.6ti;^ 


CMS [24] 


pp — > bb 




CMS [26] 


pp ZZ* i+i-i+r 




CMS [241 


pp^ww* ^i+vrv 




CMS [241 


pp — 7- 77 Untagged 


q 0+1-9 
'-'■^-1.8 


CMS [25] 


pp —7- 77 Untagged 1 


0.71?:^ 


CMS [25] 


pp — i- 77 Untagged 2 




CMS [25] 


pp — )■ 77 Untagged 3 




CMS [25] 


pp -rijj 


^■^-2.0 


CMS [25] 



Table L Processes considered in our analyses for the LHC 7 TeV run and for the Tevatron. We present the errors and best fit 
point for the signal strength for each topology. 



Channel 


exp 


Comment 


pp ->■ ZZ* £+e-£+£- 


1 q+0.6 


ATLAS @ 7 and 8 TeV [2 


pp -f-f 


1 a+0.6 


ATLAS [27| 


pp — )■ TT 


"•^-0.7 


CMS @ 7 and 8 TeV [8| 


pp — >■ bb 


n 4+0-9 

"■*-0.8 


CMS [26] 


pp ZZ* -> e+ri+r 


0.71°:! 


CMS [26| 


pp^ww* -^i+vrv 




CMS @ 7 and 8 TeV [8| 


pp — >■ 77 Untagged 




CMS [25| 


pp — i> 77 Untagged 1 


J^-0-1.0 


CMS [25] 


pp — >■ 77 Untagged 2 




CMS [25| 


pp — 77 Untagged 3 




CMS [25| 


pp — >■ 77ji loose 


U.O_2 


CMS [25| 


pp 77j j tight 


1 0+1.5 


CMS [25| 



Table IL Available data including the 8 TeV run. We present the errors and best fit point for the signal strength for each 
channel. The data that have been combined is indicated as "@ 7 and 8 TeV" . 

pulls associated with the theoretical uncertainties for the other production mechanisms as well as possible systematic 
correlated errors, however at this moment these are sub-leading effects. 

One important approximation in our analyses is that we neglect the effects associated with the distortions of 
the kinematical distributions of the final states due to the Higgs anomalous couplings arising from their non SM-like 
Lorentz structure. Thus we implicitly assume that the anomalous contributions have the same detection efficiencies as 
the SM Higgs. A full simulation of the Higgs anomalous operators taking advantage of their special kinematical features 
would increase the current sensitivity on the anomalous couplings. It would also allow for breaking degeneracies with 
those operators which only lead to an overall modification of strength of the SM vertices. But at present there is not 
enough public information to perform such analysis outside of the experimental collaborations. 

In order to predict the modification of the observables we need to include the effect of the anomalous operators in 
the production channels as well as in the decay branching ratios. As a first approximation we can assume that the K 
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Channel 


^99 






(■ttH 


pp — >■ 77 Untagged 


0.13 


0.46 


0.70 


1 


pp — ^ 77 Untagged 1 


0.57 


0.49 


0.67 


1 


pp — >■ 77 Untagged 2 


1 


0.56 


0.76 





pp — >■ 77 Untagged 3 


1 


0.56 


0.76 





PP 77J j 


0.029 


1 


0.019 






Table III: Weight of each production mechanism for the different 77 categories in the CMS analyses of the 7 TeV data. 



factor associated with higher order corrections is the same for the SM and anomalous contributions, so we write 



ano 

ay 



ay 



(9) 



where the ratio of the anomalous and SM cross sections of the subprocess Y (= gg, VBF, VH or tiH) is evaluated at 
tree level and it is multiplied by the value for the state-of-the-art SM cross section calculations (cry^^|soa) presented 
in Ref. [i^. Analogously we write the decay width into the final state X as 



\h^X) 



TS^i{h ^ X) 



T^'^'ih ^ X) 



(10) 



where the SM result T^'^\h — X)\soa is also obtained from Ref. [i^. The total width and branching ratios are 
evaluated following this recipe. We use the SM cross sections and decay widths and compute our predictions for 
niH = 125 GeV. The observed Higgs mass by ATLAS (126.5 GeV) and CMS (125.3 GeV) [3 are compatible 
within the experimental errors. We verified that the impact of changing the Higgs mass to 126.5 GeV is a sub-leading 
effect and does not alter our results. We did not include in our analyses an eventual invisible decay of the Higgs 
[13, [EH, therefore the total width is obtained by summing over the decays into the SM particles. 

The search for Higgs decaying into bb pairs takes place through Higgs production in association with a W or a Z 
so we can write 



^^bb 



'WH 



' ZH 



Br°"°[/i^6&] 



'WH 



'ZH 



Br'^^^ih^bb] 



(11) 



with the superscripts ano (SM) standing for the value of the observable considering the anomalous and SM interactions 
(pure SM contributions). 

The CMS analyses of the 7 (8) TeV data separate the 77 final into 5 (6) categories and the contribution of each 
production mechanism to a given category is presented in Table 2 of Ref. [25|. Therefore, we write the theoretical 
signal strength in these cases as 



CMS _ ^99^99 + + ^VBFCTvBF + ^VH WwH + ^ZH I + '^ttH^^ttH 



''99" gg 



eVBFCTvBF 



Br°"°[fa^77] 
Br^*^[/i^77] 



(12) 



where is the pull associated with the gluon fusion cross section uncertainties, and the branching ratio and the 
anomalous cross sections are evaluated using the prescriptions © and pUj) . The weight of the different channels to 
each category is encoded in the parameters ex with X = VBF, gg, VH, and ttH and they arc presented in Tables Hill 
andHVl 



With the exception of the above processes, all other channels are treated as inclusive, so we write the expected 
signal strength of the final state F as 



fJ-F 



99 



'(l+es)+<BF + 



'WH 



'ZH 



„ano 
^tiH 



jSM 

99 



SM 
"VBF 



^SM 
'WH 



^SM 
'ZH 



Br""°[;i^F] 



(13) 



Here we also use Eqs. © and (jlOp to obtain the anomalous cross sections and branching ratios. 

For some final states the available LHC 8 TeV data has been presented combined with the 7 TeV results. Given the 
limited available information on errors and correlations, we construct the expected theoretical signal strength as an 
average of the expected signal strengths for the center-of-mass energies of 7 and 8 TeV. We weight the contributions 
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Channel 


^33 


tVBF 




(■ttH 


pp — >■ 77 Untagged 


0.11 


0.25 


0.48 


1 


pp — ^ 77 Untagged 1 


0.59 


0.50 


0.72 


1 


pp — >■ 77 Untagged 2 


1 


0.54 


0.58 





pp — >■ 77 Untagged 3 


1 


0.54 


0.78 





pp — i> 77jj loose 


0.094 


1 


0.064 





pp 77j j tight 


0.024 


1 









Table IV: Same as Table |III] but for the 8 TeV CMS data. 



by the total number of events expected at each energy in the framework of the SM, i.e. given a final state X we 
evaluate 

,7TeV „SM,7Tc'V p , ,,8TeV ^SAf,8TeV , 



comb „ '^X J-rTcV + ^J■X X ^8TeV , . 

SM,7TcV r , SA/,8TcV ' ^ ) 

-i-YTcV + ""jf i-8TcV 



where 'C7(8)tgV stands for the integrated luminosity at 7 (8) TeV accumulated in the channel being analyzed. When 
considering the full available data set we consider all the processes in Table U and [ITl neglecting the LHC 7 TeV 
processes whose data has been combined with the 8 TeV run; we indicate in Table |TT] that the data has been combined 
by "@ 7 and 8 TeV". 

The evaluation of the relevant tree level cross sections was done using the package MadGraphS [s^l with the 
anomalous Higgs interactions introduced using FeynRules [s^. We also cross checked our results using COMPHEP [HI, 
Issj l and VBFNLO [sBl. The evaluation of the partial width was done using the expressions presented in Ref. [s^ . 



IV. CONSTRAINTS ON THE HIGGS ANOMALOUS INTERACTIONS 

We next derive the allowed values of the Higgs interactions to vector bosons using the available Tevatron data [l^l , 
ATLAS 7 TeV and 8 TeV 0,[23 results, and CMS 7 TeV [Hi^ and 8 TeV [i,[2l,[2^ data. The results are 

presented in Fig. [T]-Fig. [5] where several one-dimensional and two-dimensional projections of the Ax^ = _ 
function(s) are shown. We find Xmm = 12.12 (12.13) in scenario I (II) for the global analysis (i.e. for a total number of 
26 data points). The SM lays at Xsa/ = 20.87 i.e. within the 96.7% (98.7%) CL region in the three- (two-) dimensional 
parameter space. The corresponding best fit values and 95% allowed ranges are summarized in Table IVl 

Figure [1] shows Ax^ as a function of each of the operator coefficients in scenario I after marginalizing with respect to 
the two undisplayed ones. To illustrate the effect of the different data sets, the results are shown for three combinations 
of the available data: the dotted (dashed) line stands for the results obtained using only the LHC 7 TeV (LHC 7 
TeV and Tevatron) data while the solid line is derived using the full available data set. The central panel of this 
figure displays the Ax^ dependence on j^. As we can see the analysis of the LHC 7 TeV data only leads to a large 
flat region around the minimum indicating that this data set has a small sensitivity to /vi/, «-e. the Higgs couplings 
to W and Z pairs. This is expected since the 77 channel is the dominant observable in this sample. The addition 
of the Tevatron data, dominated by the Higgs associated production, enhances the sensitivity to deviations in HZZ 
and HW^W^ , that is, to smaller values of fw- The addition of the LHC 8 TeV results further tightens the allowed 
values giving, for the global analysis, the constraint — 13 < fw < 20 at 95% CL. 

Ax^ as a function of fg is shown in the left panel of Figure [T] where we see that the analysis present two totally 
degenerate minima leading to two distinct allowed ranges. This degeneracy - as others that we encounter in this work 
- is due to the interference between the SM and anomalous contributions. We see that before the inclusion of LHC 8 
TeV data, the two allowed ranges overlapped at CL higher than 90 %, while in the global analysis they are separated 
at more than 3a. The value of the gluon fusion cross section at the minima is around 43% of its SM value (see left 
panel of Fig. U]) and this cross section is highly suppressed in the region between the minima. This highly suppressed 
gluon fusion cross section was not completely disfavored before the 8 TeV data because the CDF & D0 and CMS @ 
7 TeV bb channels still allowed for a larger fw coupling to enhance associated production, compensating the large 
reduction of the gluon fusion cross section (see also discussion of Fig. [3]). After the inclusion of the LHC 8 TeV data 
this is no longer possible. So altogether the global analysis constrains fg to lie in one of the two intervals [—0.3 , 7.3] 
or [15 , 23] at 95% CL. 

The Ax^ dependence on fww in the scenario I is presented in the right panel of Fig. [TJ A salient feature of this 
plot is that Ax^ is concentrated around two narrow non-overlapping regions centered around almost (but not totally) 
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degenerate minima. Unlike for fg, the two minima in fww = fsB are not fully degenerated because these operators 
modify not only the Higgs coupling to photons but also to WW and ZZ and the contributions to these last two 
vertices are slightly different at the two minima. Moreover, we also see that Tevatron data has a limited impact on 
this parameter while the inclusion of the LHC 8 TeV results tighten the bounds on fww which at 95% CL is bounded 
to lie in one of the two intervals [—0.8 , —0.1] or [1.5 , 2.2]. 




Figure 1: The left (central, right) panel exhibits Ax^ as a function of fg [fw, fww) in the framework of scenario I. Each panel 
contains three lines: the dotted (dashed) line was obtained using only the LHC 7 TeV (LHC 7 TeV and Tevatron) data while 
the solid line stands for the result using all the available data. In each panel Ax^ is marginalized over the two undisplayed 
parameters. 

The dependence on the scenario considered is illustrated in Figure [2] where we plot the Ax^ dependence on fg and 
fww of the global analysis in scenarios I and II. As we can see the results for both scenarios are almost coincident in 
both panels. This is due to the fact that in scenario I the full available data set is well described by fw = /b — for 
all allowed values of fww and /g, consequently the two scenarios give very similar results. 

Let us turn our attention towards the correlations between the three free parameters of scenario I. Figure [3] depicts 
68%, 90%, 95%, and 99% CL (2dof) allowed regions of the fww ® fg (upper right panel), fw ® fg (upper left panel) 
and fw fww (lower panel) planes using all attainable data. We obtained these plots marginalizing over the free 
parameter not appearing in each of the panels. 

We can see from the upper right panel of Fig. [3] that there are four well isolated allowed "islands" in the fww fg 
plane. Moreover, within each of these islands fww and fg are strongly correlated or anti-correlated. As mentioned 
before, the existence of degenerate islands is due to the interference between the SM and anomalous contributions 
which allow two different values of the anomalous couplings to lead to the same cross section or branching ratio. In 
the case at hand, the gluon fusion cross section preferred by the fit is around 43% of its SM value. It is interesting to 
notice that if the results from the bb channel are removed from the fit the vertical gap between the two islands on the 
left (or on the right) disappears - that is, intermediate values of fg (which correspond to further suppressed gluon 
fusion production) become allowed. This happens because the bb data, which is dominated by associated production, 
constrains the coupling of the Higgs to W and Z pairs. In our framework, this leads to (a) an associated upper bound 
on the Hjj branching ratio, and (b) an upper bound on VBF and associated production. 77 data mainly restricts 
the product of the gluon fusion cross section and the Higgs branching ratio into photons, thus weakening the upper 
bound on the latter allows the former to have smaller values. Furthermore even smaller gluon fusion cross sections 
are permited because of the possible increase in the VBF and associated production processes. 

The upper left panel of Figure [3] shows the presence of two isolated regions in the fw (8> fa plane and that there 
is a very weak correlation between the parameters within each region. Here again, the removal of the bb data leads 
to the disappearance of the gap between the allowed regions. The lower panel displays a behavior similar to the one 
observed in the upper left, but in the fw €5 fww plane. 

For the sake of completeness we also show the results of the global analysis in scenario I in terms of the allowed 
ranges of Higgs production cross sections and decay branching ratios in Fig. 2] and Fig. [51 The results shown in these 
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Figure 2: Ax^ as a function of /g (left panel) and fww = fsB (right panel) for the full combined analysis. The solid lines 
correspond to scenario I in which Ax^ is marginalized over the two undisplayed parameters in each panel: fww = /ss and 
fw = /s in the left panel, and fg and fw = /s in the right panel. The dashed lines correspond to scenario II, i.e. imposing 
first the prior fw = fs ~ and then marginalizing over fww = fsB (left) and fg (right). 

figures are obtained by projecting the three-dimensional A^;^ function on tlic displayed observables and marginalizing 
on the independent undisplayed combination(s). 

Finally we also verified that the results do not change significantly when we do not employ the pulls to perform the 
fit. This behavior could be anticipated since the experimental errors are still much larger than the errors described 
by the pulls; a situation that will change as more statistics accumulate. 



V. DISCUSSION 

Once a Higgs boson like state has been discovered we must study its properties to establish if it is indeed the 
state predicted by the SM. In addition to that, it is also important to look for additional states that might play 
a role in the electroweak symmetry breaking. In this article we have studied the Higgs couplings to gauge bosons 
using a model-independent characterization of the deviations with respect to the SM values in terms of dimension-six 
operators and the available data from Tevatron and LHC at 7 TeV and 8 TeV. This approach still assumes that the 
Higgs field is a doublet of the SU{2)l symmetry and that the deviations of its couplings from the SM values are due 
to additional heavy states. Notwithstanding, our framework allows for independent modifications of the couplings to 
gluons, photons and weak gauge bosons. 

In this study we have demonstrated that the present available data is enough to start gaining some information on 
the different Higgs couplings to gauge bosons. For instance, our analyses indicate that a reduced gluon fusion cross 
section is preferred when we use the full available data set, with the most favored value being 43% of the SM value. 
We can see this preference for a reduced gluon fusion cross section in the right panel of Figure |4] while the VBF and 
associated production cross sections are in agreement with the SM prediction. From this panel we can extract that 
the 95% CL allowed region of the gluon fusion cross section is [0.1, 1.1] times the corresponding SM value. This is 
consistent with the CMS analyses [8| which, using a different framework, also points in this direction as a reduced 
coupling of the Higgs to top quarks is preferred by their results. 

Taking into account that the presently measured 77 yield is above the SM prediction, the diminished gluon fusion 
cross section points to an enhanced Higgs branching ratio in 77; a fact that can be observed in our analyses. The left 
panel of Figure H] shows that the 77 branching ratio is indeed augmented, with a best fit value of 2.9 times the SM 
value and the 95% CL allowed region being [1.4,5.4] times the SM branching ratio. Furthermore, we can see from 
this panel that the Higgs branching ratio into W~^W~ and ZZ is in agreement with the SM expectations. 

Presently the 77 channel is the best measured channel and its rate is above the SM prediction. The operators 
Oww, Obb and Ogg are the ones affecting this channel, therefore they are the ones showing the largest impact of 
the full data set. This can be seen from the strong correlations and the well isolated islands present in the upper right 
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Figure 3: 68%, 90%, 95%, and 99% CL (2dof) allowed regions of the plane fww fg (upper right panel), fw ® fg (upper left 
panel) and fw ® fww (lower panel) using all available data. These results are obtained for scenario I and after marginalization 
over the undisplayed parameter in each panel. The best fit points are indicated by a star while the second local minima are 
indicated with a dot. 

panel of Figure [3] as well as by the correlations between the gluon fusion cross section and the Higgs branching ratios 
into electroweak gauge bosons in Figure [5] From the upper left panel of this figure we can see an anti-correlation 
between the gluon fusion cross section and the Higgs branching ratio into two photons; once again it is clear that 
there is a preference for reduced gluon fusion cross sections and enhanced decay into photon pairs. The other two 
panels of Figure [S] show the mild dependence of the Higgs branching ratio into W~^W~ with the gluon fusion cross 
section or the two photon branching ratio. 

Our analyses of scenario I also shows that the presently available data prefers small values of fw = Jb', see the 
central panel of Figure [TJ This indicates that large deviations in HZZ and HW~^W~ interactions, as well as to triple 
gauge-boson couplings, are not favoured. This behavior was expected because the data points for Higgs couplings 
to W^'s and Z's are in agreement with the SM within la; see Figure S] left panel. Furthermore, the present direct 
constraints on triple gauge-boson vertices lead to bounds on fw that are of the same order as the ones derived 
here from Higgs phenomenology. So in the future the combined analysis of Higgs data and measurements of the 
anomalous triple gauge-boson couplings can be used to reduce the degeneracies observed in our results since they 
present a different dependence on the anomalous couplings fw and /b; see Eqs. (j?]) and ([H]). In this respect, it 
is interesting to notice that electroweak precision measurements still give rise to the tightest limits on the Higgs 
anomalous interactions [sO, [111 . 

We finish with a word of warning. The precise numerical results presented here, that arc summarized in Tabic IVl 
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Figure 4: Ax^ as a function of Higgs branching ratios into electroweak gauge bosons (left panel) and the cross section for 
different production processes (right panel) normalized to the SM values. In the left panel the solid (dashed, dotted) line stands 
for the branching ratio into 77 {W^W~ , ZZ), while, in the right panel, the solid (dashed, dotted) line represents the gluon 
fusion (VBF, VH) production cross section. 



LHC (7+8)TeV + Tevatron 




DrS^ / DKy^ 

Figure 5: Allowed regions for several combinations of Higgs branching ratios and production cross section. In each panel Ax^ 
is marginalized with respect to the combination of couplings independent of the two displayed observables. As in Figure [3] the 
regions are shown at 68%, 90%, 95%, and 99% CL (2dof). 



should be taken with a grain of salt; due to the simplifying hypothesis used in our analyses we should be aware that 
details can change if a more complete approach is used. Nevertheless we verified that our results arc rather robust 
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Best fit 


95% CL allowed range 


fw = /s (TeV-2) 


-0.8 


[-13,20] 


fww = fsB (TeV-2) 


-0.4, (1.8) 


[-0.8,-0.1] and [1.5,2.2] 


(TeV-2) 


3.7, 19 


[-0.3, 7.3] and [15, 23] 




2.9 


[1.4,5.4] 


D ryano i jj rySM 


1.1 


[0.8, 1.3] 


BR%"£/BR§%' 


1.1 


[0.7, 1.3] 


ano 1 SM 
99 1 99 


0.4 


[0.1,1.1] 


ano 1 SAI 
^VBF / ^VBF 


1.0 


[0.8, 1.5] 


ano 1 SlS.1 


1.1 


[0.6,2.1] 



Table V: Best fit values and 95% CL allowed ranges for the combination of all available data. For /g we show the two degenerate 
best fit values. For fww = /bb together with the best fit we show in parenthesis the value at the second minimum. 



when we use only parts of the available data. 
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